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We compute the K-edge indirect resonant inelastic x-ray scattering (RIXS) spectrum of a triangular lattice 
antiferromagnet in its ordered coplanar 3- sublattice 120° magnetic state. Conventional K-edge RIXS spectra 
prohibits the presence of odd spin flip terms. However noncollinearity of the spin arrangement in a triangu¬ 
lar lattice causes the transverse and longitudinal spin components to be coupled giving rise to intrinsic odd 
spin flip trimagnon excitations. By considering the first order self-energy corrections to the spin wave spec¬ 
trum, magnon decay rate, bimagnon interactions within the ladder approximation Bethe-Salpeter scheme, and 
the effect of three-magnon contributions up to 1/S - order we find that the RIXS spectra is non-trivially af¬ 
fected. For a purely isotropic triangular lattice model, the peak splitting mechanism and the appearance of a 
multipeak RIXS structure is primarily dictated by the damping of magnon modes. At a scattering wavevector 
corresponding to the zone center T point and at the roton point q = M where the magnon decay rate is zero 
a stable single peak forms. However, the microscopic origins of these peaks are different. At the T point, the 
contribution is purely trimagnon at the 1/S level and occurs approximately at the trimagnon energy of 6 JS. 
This provides experimentalists with a means to detect purely trimagnon excitations, even at the K-edge. The 
roton peak occurs at a lower energy of 4 JS. The K-edge single peak RIXS spectra at the roton momentum 
can be utilized as an experimental signature to detect the presence of roton excitations. A unique feature of 
the triangular lattice K-edge RIXS spectra is the nonvanishing RIXS intensity at both the zone center T point 
and the antiferromagnetic wavevector K point. This result is in sharp contrast to the vanishing K-edge RIXS 
intensity of the collinear ordered magnetic phases on the square lattice. We find that including XXZ anisotropy 
leads to additional peak splitting, including at the roton scattering wavevector where the single peak destabilizes 
towards a two-peak structure. The observed splitting is consistent with our earlier theoretical prediction of the 
effects of spatial anisotropy on the RIXS spectra of a frustrated quantum magnet [Luo, Datta, and Yao, Phys. 
Rev. B 89, 165103 (2014)]. In summary, the features of an indirect K-edge RIXS spectra of a triangular lattice 
quantum magnet can be interpreted as a combination of magnon decay and spin anisotropy effects. 

PACS number(s): 78.70.Ck, 75.25.-J, 75.10.Jm 


I. INTRODUCTION 

Contrary to the historical prediction of the spin- ^ tri¬ 
angular lattice antiferromagnet (TLAF) as a canonical ex¬ 
ample of a spin liquid state [1], extensive - theoretical [2- 
4], numerical [5-11], and experimental [12-20] - studies on 
the nearest-neighbor Heisenberg model has established the 
ground state configuration as a 120° long-range coplanar 3— 
sublattice arrangement. The predicted ordering pattern per¬ 
sists for all values of spin S , including the S — \ state 
where quantum fluctuations lead to a 60% suppression of the 
magnetic order parameter from its classical Neel value [5- 
8], At present there exists a plethora of real TLAF materials, 
with both isotropic and anisotropic interactions, which pro¬ 
vide a motivation to study triangular lattice frustrated mag¬ 
nets [12, 15-21], Further impetus to investigate and de¬ 
lineate the physical properties of the TLAF stems from the 
flurry of recent theoretical and numerical investigation to clar¬ 
ify the ground and excited state properties of both isotropic 
and anisotropic triangular lattice systems [22-31], 

Traditionally, information on the magnetic ground state 
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and single-magnon excitations is inferred from inelastic 
neutron scattering (INS) experiments [32, 33], However, 
with enhancements in instrumental resolution of the next 
generation synchrotron radiation sources resonant inelastic 
X-ray scattering (RIXS) spectroscopy offers the condensed 
matter and materials science community an alternate op¬ 
tion to experimentally probe magntic excitations in corre¬ 
lated magnets [34], As a photon-in photon-out spectro¬ 
scopic technique, RIXS can offer direct information on both 
single-magnon and multimagnon excitations. Present efforts 
to understand the K- edge indirect RIXS spectra are primar¬ 
ily directed towards the study of square lattice compounds 
in the Neel antiferromagnetic and collinear antiferromagnetic 
phases [35-41]. In a recent publication. Ref. [40], the authors 
of this paper have shown that in the case of an anisotropic 
square lattice with strong frustrating further neighbor inter¬ 
actions the RIXS spectrum can split into a robust two-peak 
structure, over a wide range of transferred momenta, in both 
magnetically ordered phases. It was also predicted that the 
unfrustrated model contains a single-peak structure. 

In RIXS spectroscopy single and three spin flip processes 
are allowed at the L- and M- edges due to the presence of 
spin-orbit coupling [42-44], But, in a square lattice, excita¬ 
tions of odd magnons are prohibited at the K-edge [34] and 
the spectra originates purely from the bimagnon contribution. 
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In the absence of an external magnetic field the spin order¬ 
ing in a square lattice system is collinear and the magnon 
excitations are long lived without any damping. In contrast, 
in the TLAF the noncollinear ground state contains inherent 
three-magnon excitations (odd spin flip terms). The coupling 
of the longitudinal and transverse spin excitations gives rise 
to a finite quasi-particle lifetime (see Fig. 1) which intro¬ 
duces an intrinsic damping of the magnon modes [45, 46]. 
Hence, the presence of the trimagnon interaction, even at the 
K-edge, motivates several unanswered questions within the 
context of quantum magnetism and RIXS spectroscopy. How 
does the presence of an intrinsic damping affect the indirect 
K-edge RIXS spectra? What role does the interplay between 
geometrical frustration and spin anisotropy have on the RIXS 
spectra? In this article, we predict the effects of bimagnon and 
trimagnon processes in indirect RIXS spectroscopy of a geo¬ 
metrically frustrated TLAF, a topic which is unexplored both 
theoretically and experimentally. 

The microscopic mechanism underlying magnetic excita¬ 
tions in the indirect RIXS process involves a local modifica¬ 
tion of the superexchange interaction mediated via the core 
hole [47, 48]. The resulting RIXS spectra is expressed as a 
momentum-dependent four-spin correlation function which 
can probe bimagnon excitations across the entire Brillouin 
zone (BZ) [37-39]. Hence, RIXS is complementary to optical 
Raman scattering which is restricted to zero momentum [49- 
52]. From a theoretical perspective elucidating the nature of 
the bimagnon dispersion affected both by two-magnon ladder 
scattering processes and three-magnon interactions is chal¬ 
lenged by the appearance of several non-trivial dynamical 
properties in the spin-wave excitation spectrum [46, 53-55]. 
Namely - (i) strong renormalization of magnon energies with 
respect to the linear spin-wave theory result, (ii) finite life¬ 
time due to spontaneous magnon decays at zero temperature, 
and (iii) appearance of rotonlike minima at the edge center of 
the BZ (M point, see Fig. 1(a)). 

The objective of this paper is to elucidate the role of 
magnon-magnon interaction, spontaneous magnon decays, 
the effect of the rotonlike minima, and spin anisotropy on 
the indirect RIXS spectra of a TLAF. For this purpose, we 
consider an easy-plane XXZ triangular lattice model. In the 
isotropic limit, the XXZ model can provide an accurate de¬ 
scription of the Ba 3 CoSb 2 C >9 system [12, 56]. In addition, it 
provides a starting point for the discussion of RIXS effects in 
anisotropic TLAF [15-20], We compute the RIXS intensity 
utilizing the 1/5 spin-wave expansion technique within the 
Bethe-Sapleter scheme where interaction effects arising from 
both the quartic terms via the ladder scattering process and the 
contributions of the cubic anharmonic terms up to 1/5 order 
are included. 

The main results of our article can be summarized as fol¬ 
lows. First, in the case of an isotropic nearest neighbor TLAF 
we find that the spontaneous magnon decay and kinematic 
constraints of the phase space inherent to the model is the 
primary cause for creating a multipeak (more than two-peak) 
structure in the RIXS spectra. Second , contrary to the K-edge 
RIXS intensity of the square lattice case, in the TLAF the 
RIXS intensity does not vanish at the I point and at the K 


point. At the F point, the bimagnon intensity is zero and the 
single peak spectra results purely from the trimagnon con¬ 
tribution, approximately at energy scale of 6JS correspond¬ 
ing to the three magnon energy. This provides experimental¬ 
ists with a means to detect purely trimagnon RIXS spectra at 
the K-edge. At the antiferromagnetic wave vector K point 
the RIXS intensity is dominated by the bimagnon excitations. 
Third, an important conclusion of our work is the proposal of 
utilizing RIXS as a probe to detect the presence of the roton 
mode. We show that at a scattering wave vector equal to the 
roton momentum q — M the RIXS spectra has a single peak 
structure. Barring the T point peak which occurs at a higher 
energy, at all other special high symmetry points of the mag¬ 
netic BZ the RIXS spectra splits into a multipeak structure. 
The appearance of the single peak structure can serve as an ex¬ 
perimental signature to detect the appearance of a roton mode 
in a TLAF. Fourth, including the XXZ anisotropy leads to fur¬ 
ther peak splitting including at the roton scattering wavevec- 
tor point. Fifth, we find that the conceptual signature of slow 
moving bimagnons as an indicator of RIXS peak splitting (in¬ 
stability), as proposed in our earlier work on the two-peak 
splitting theory within the context of the anisotropic square 
lattice Heisenberg model [40], still holds (see Fig. 8). 

This article is organized as follows. In Sec. II, we introduce 
the XXZ Hamiltonian, present the expression for the effective 
Hamiltonian within an interacting spin-wave formalism, and 
compute the intensity maps for the renormalized energy and 
magnon decay rate, up to 1 /S corrections. In Sec. Ill, we state 
the definition and the expression of the TLAF RIXS scattering 
operator containing both the bimagnon and trimagnon contri¬ 
butions. In Sec. IV, we display our results, state the formal¬ 
ism and numerical approach for computing RIXS intensity, 
and discuss the implications of our result within the context 
of a TLAF (geometric frustration). First, in Sec. IV A, we 
present the results for the noninteracting bimagnon and tri¬ 
magnon RIXS intensity and spectral weight. In Sec. IV B, we 
outline our formalism and calculate the interacting bimagnon 
intensity. In Sec. IV C, we calculate the full RIXS spectrum. 
In Sec. V, we present our concluding remarks and discuss the 
appearance of slow moving bimagnons as a signature of peak 
splitting. Finally, to preserve clarity in the main body of the 
text, we state the details of the spin-wave theory derivation of 
the effective Hamiltonian in Appendix A and display results 
to validate our numerical approach in Appendix B. 


H. TLAF MODEL AND MAGNON DECAY 

Inelastic neutron scattering data of a TLAF reveals 
well-defined sharp modes in the low-energy excitation 
spectrum accompanied with a broad continuum at inter¬ 
mediate and high energies [56-58]. A number of com¬ 
peting theoretical proposals, ranging from a proximate 
spin-liquid phase [25, 28] to enhanced magnon-magnon in¬ 
teractions [59-61] have been proposed to explain the nature 
of the spin wave excitation spectrum. Our starting point is 
the spin S, nearest-neighbor XXZ antiferromagnetic model 
on the triangular lattice. The spin-wave theory Hamiltonian 
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FIG. 1. (Color online) Intensity maps of the 1/S spin wave spectrum 
for the S = 1/2 triangular-lattice antiferromagnet with easy-plane 
anisotropy a = 1 (upper row) and a = 0.95 (lower row), (a) and (c), 
renormalized magnon energy e k . (b) and (d), magnon decay rate r k . 
r = (0,0), E = (2n/3, 0), X = Or, 0), K = (4zr/3, 0), M = (0,2 n/ V3), 
and Y = (0, n) V3) points in (a) are highlighted. Note, the over¬ 
all damping of spin waves is strongly reduced and the decay region 
shrinks with increasing anisotropy. The realistic magnon decays dis¬ 
appear at a x 0.92 [55], 


in the local (x - z ) rotating frame associated with the ordering 
wave vector Q = ( 4n/3,0 ), K— point in BZ, takes the follow¬ 
ing form [55]: 

<H =jY j [aS y i S y j + cos (6i - 0,)(5j5) + 5*5*) 

(O'> 

+ S m(e i -9j)(S z i S x j -S*S z j )]. (1) 

where 6, = Q ■ r, and we have also introduced an easy-plane 
anisotropy parameter a e [0,1]. In Appendix A we outline the 
derivation of the effective interacting spin-wave Hamiltonian 
'He sf in the first order 1/5 expansion with respect to linear 
spin-wave theory. The resulting expression is 

H e ff =^Ok + 5s k )b[}j k + ^ ^ V a (b\b\bi + H.c.) 

k |k,| 

+ i Z V b (b \ b l b l + H - c -) + Z V ' b l b l b 3^4- (2) 

Ik,I {k,t 

where we have adopted the convention that 1 = ki, 2 = k 2 , 
etc, and momentum conservation is assumed for various k- 
summations. The bare magnon dispersion given by the linear 
spin-wave theory is expressed as 

e k = 3 JS V(l-n)d+2Q7k), (3) 

with y k = |[cosk f + 2cos(k x /2)cos(y/3ky/2)]. The explicit 
forms for the interacting vertices 6s, V a j,, and V c are detailed 
in Appendix A. At zero temperature the bare magnon propa¬ 
gator is defined as 

Gg'fk, oj) = oj- s k + ;0 + . (4) 


The first order 1/5 correction to the magnon energy is deter¬ 
mined by the Dyson equation 


<jj - s k - Z(k, oj) = 0, 


(5) 


with the one-loop self-erengy X(k, of) = X„(k, to) + £/,( k, to) + 
Z f (k), where X r (k) = is a frequency-independent Hartree- 
Fock correction, while &(k, oj) are calculated as [4, 46, 54, 
55] 


E fl (k,cu) = ^ 

P 


|y f ,(p,k-p;k)| 2 
to - s p - e k _ p + ;0 + ’ 


S i (k,oj) = -^, 

P 


|y fa (p, -k - p, k)[ 2 

to + e p + e k+p - ;0 + 


( 6 ) 

(7) 


The on-shell solution consists of setting to - s k in the self¬ 
energy (6) and (7) leads to the following expression for the 
1/5 renormalized spectrum 


w k = w k - /F k - e k + E(k, s k ). (8) 

In Fig. 1 we display the intensity maps for the renormaliza¬ 
tion of magnon energy e k and the magnon decay rate F k for 
the 5 = 1/2 triangular antiferromagnet for a - 1 (isotropic) 
and a = 0.95 (anisotropic) case. From Fig. 1(d) we observe 
that the magnon decay rate decreases drastically in the pres¬ 
ence of anisotropy. This is due to the reduced phase volume 
where the kinematic constraint e k = s p - e k _ p in the self¬ 
energy (6) is satisfied. The magnon decay intensity maps 
Fig. 1(b) and Fig. 1(d) play an important role in understanding 
the origins of the multipeak RIXS structure shown in Fig. 6 
and Fig. 7. 


HI. INDIRECT RIXS CORRELATOR 

In Mott insulating systems, multimagnon excitations can be 
created dynamically by the presence of the core-hole potential 
in the intermediate state of indirect RIXS process. The effec¬ 
tive scattering operator, in the first order, under the assumption 
of the ultra-short core-hole life-time (UCL) expansion is given 
by [37, 39] 


<R q = jY J e iqr ‘SrS i+s , (9) 

i ,6 

where r, is the position of the ion absorbing the incident pho¬ 
ton and 6 denotes the neighboring vectors. After consecu¬ 
tive Holstein-Primakoff and Bogoliubov transformations, the 
magnon creation parts of the RIXS scattering operator can be 
expressed in terms of the bosonic operators as 

<R q = Yj M(l,2)b\bl+ Yj N(\,2,2>)b\b\b\, ( 10 ) 

1+2=q l+2+3=q 

with the bimagnon scattering matrix element expression is 
given by 

3/5 , 1 

4/(1,2) +r q + (0'- + y2)](u\V2 + V\U 2 ) 
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-(<* + ^Xyi +yi)(uiu 2 + U1U2)}, (li) 

and the trimagnon scattering matrix element is given by 

3 js I 3 1 

N(l,2,3) = — iyj — [(Ti -72+3 - ^rqKwi + «l) 

X(m 2 U3 + V 2 U 3 ) + (j 2 ~ 71+3 - ^7q)(«2 + V 2 ) 

X(mv 3 + V 1 U 3 ) + (73 - yi +2 - ^7q)(“3 + V 3 ) 

x(uiv 2 + viu 2 )], (12) 

where u, v and 7 are defined in Appendix A. The three-boson 
term in our theory has no analog in the collinear phases of 
a square lattice quantum magnet. Note, the corrections from 
magnon interactions for the trimagnon intensity appear at the 
1/5 2 order and are neglected in the remainder of this paper. 

The frequency- and momentum-dependent magnetic scat¬ 
tering intensity is related to multimagnon response function 
via the fluctuation-dissipation theorem. The full 1/5 cor¬ 
rections to the indirect RIXS susceptibility is the sum of bi- 
magnon and trimagnon contributions given by 

/(q,w) = --Im[y 2 (q,w) +* 3 (q,w)], 

71 

- --Im[y RIX s(q.w)], (13) 

n 

which involve an interacting two-magnon susceptibility 
^f 2 (q, oj) and a non-interacting three-magnon susceptibility 
;^(q, oj). The susceptibilities can be expressed explicitly from 
the corresponding multi-magnon Green’s function defined as 

* 2 (q, w) = 2 M k M k 'Yl kk (q, oj), (14) 

kk' 

X3 (Q? ^) = ^ ] ^k,p^k',p'^kp;k'p'(q? (15) 

kp;k'p' 

where If and A are denoted as the bimagnon and trimagnon 
propagators, respectively. The momentum-dependent two- 
magnon and three-magnon Green’s function in terms of Bo- 
goliubov quasiparticles are defined as 

iUkk-ta, t) = (Tb kl<l (t)b kWv.q/7 k ,), (16) 

«A k p; k 'p'(q, t) = <T'/ 7 k ( 0 Vk-p( 0 ( 2 p( 0 ^( 2 q_ k ,_p^p-),( 17 ) 

where T is the time-ordering operator and (•) is the average 
of the ground state. In the following sections, using Eq. (16) 
and Eq. (17), we will compute the noninteracting and the in¬ 
teracting RIXS spectra. 


IV. RESULTS AND DISCUSSION 

A. Noninteracting bi- and trimagnon spectra 

Using Eqs. (11)— ( 12) and applying Wick’s theorem to 
Eq. (16) and Eq. (17), we obtain the following expressions 



u/3JS u/3JS u/3JS 


FIG. 2. (Color online) Non-interacting bimagnon (/ 2 (q, w)) and tri¬ 
magnon (/3 (q, ai)) RIXS intensity of the isotropic TLAF at momen¬ 
tum transfer wavevector q corresponding to the special high symme¬ 
try points of a triangular lattice magnetic BZ. The bimagnon DOS 
D 2 ( q, oj) is also shown as a dashed line. The nonzero intensity at the 
T and at the K point is a unique RIXS feature of the noncollinear 
ground state configuration. 



FIG. 3. (Color online) The non-interacting bimagnon and trimagnon 
total spectral weight for the anisotropic TLAF with (a) a = 1 and (b ) 
a = 0.95. Irrespective of the presence of anisotropy, the bimagnon 
and trimagnon intensity complement each other at the zone center 
and at the zone boundary. Introduction of anisotropy removes the 
singularity due to the opening up of a gap in the spin wave spectrum. 


for the noninteracting bimagnon (/ 2 (q, oj)) and trimagnon 
(/ 3 (q, a>)) scattering intensity 

h( q, w) = 2 ^ M^ +q _ k 6(oj - s k+q - £ k ), (18) 

k 

h( q, w) = 6 'Yj iV k,q-k-p.p^ W - fik - e q-k-p - Sp). (19) 

k.p 

In Fig. 2, we show the results for the 5 = 1/2 isotropic 
Heisenberg model at various points in the BZ. At the T 
point the contribution is purely from the trimagnon excita¬ 
tions, see Fig. 2(a). The bimagnon RIXS intensity displays 
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(a) (b) 




FIG. 4. (Color online) Pure trimagnon contribution to the RIXS in¬ 
tensity. (a) no anisotropy (a=l) (b) with anisotropy (a=0.95). 


a nonzero elastic peak at the K point, see Fig. 2(f). The 
indirect RIXS spectra even at the noninteracting level, in a 
noncollinear quantum magnet, exhibits significant differences 
from the collinear ordered quantum magnets where the inten¬ 
sity vanishes at the BZ center and at the antiferromagnetic 
wavevector[37-40]. 

The noninteracting bimagnon RIXS intensity in Eq. 18 is 
proportional to the bare two-magnon density of states (DOS) 

D 2 ( q, 0)) = ^ 8(a> - s k+q - s k ), (20) 

k 

A close inspection on the DOS in Fig. 2 shows that these Van- 
Hove singularities which originate from the maximum or sad¬ 
dle points of the two-magnon continuum et+q + £k partially 
transfer to the RIXS intensity, see Fig.2(b-e), and the spec¬ 
trum line shape at q = M (Fig.2(e)) resembles the DOS with 
minimal RIXS matrix element effects. 

In Fig. 3 we show the variation of the total spectral weight 
across the BZ for the bimagnon and trimagnon component, re¬ 
spectively. By using the bare intensity Eqs. (18) and Eqs. (19) 
we obtain 

/-»oo 

*S 2 (q) = I dcul 2 (q, oS) = 2 V M 2 (k + q, -k), (21) 

Jo V 

/"»00 

<S 3 (q) = I dw/ 3 (q, oS) = 6 V /V 2 (k, q - k - p; p), (22) 
Jo V 

In general, the trimagnon excitation dominates the indirect 
RIXS total spectral weight in the vicinity of the BZ cen¬ 
ter, while the bimagnon spectral weight becomes overwhelm¬ 
ingly large at the boundary of the BZ where the three-magnon 
intensity is negligible. The most remarkable feature of the 
isotropic model, see Fig. 3(a), is the elastic peak at the anti¬ 
ferromagnetic wave vector which resembles the longitudinal 
dynamic structure factor probed by neutron-scattering exper¬ 
iments [61-63], Upon inclusion of anisotropy, a = 0.95, the 
elastic peak at q = K disappears, see Fig. 3(b), since a gap is 
now introduced in the spin-wave dispersion (4) at the order¬ 
ing wave vector. 

In Fig. 4 we show the pure trimagnon contribution along the 
T —> M path obtained using the noninteracting three-magnon 
susceptibility ^ 3 (q, a>) in Eq. (15). We plot the spectra both 



FIG. 5. (Color online) Diagrammatic representation for (a) two- 
magnon propagator Il kk -(q, a>), (b) Bethe-Salpeter euqation of the 
vertex function r kk /(w, a>') and (c) the 1/S order irreducible inter¬ 
action “Vir. Solid lines with an arrow in (a) and (b) stand for the 
single-magnon propagators. 


in the presence and in the absence of anisotropy. We observe 
that the trimagnon spectra peak occurs at a higher energy ap¬ 
proximately around 6 JS around the F point, which downshifts 
before undergoing an upward shift to 6 JS around the M point. 
In the presence of anisotropy, see Fig. 4(b), there is an overall 
upward shift of the energy peak. The observed effect could 
be an artifact of considering a noninteracting trimagnon spec¬ 
tra. In the next section we consider the interacting bimagnon 
RIXS intensity up to 1/S order. 


B. Bimagnon excitations: 1/5 corrections 


We now proceed with the analysis of 1 /S correction to the 
two-magnon Green’s function by taking into account both the 
self-energy correction to the single magnon propagator G ac¬ 
cording to the Dyson equation and the vertex insertions to the 
two-magnon propagator n which satisfies the Bethe-Salpeter 
(BS) equation [64, 65]. The diagrammatic representation of 
such procedure is depicted in Fig. 5(a) and 5(b). The to¬ 
tal irreducible bimagnon scattering vertices in Fig. 5(c) fall 
into two categories, which we term as direct ('V 4 ) and indi¬ 
rect CV a f d ). 

The direct collision between the two main magnons is 
caused by the quartic vertex 'V 4 while the cubic vertices 
'V“ d represent the indirect magnon-magnon interactions. 
Note that in the direct ladder interaction events the two main 
magnons created in the RIXS process are stable while virtual 
decays and recombination are allowed in the indirect colli¬ 
sion process. Using Feynman rules in momentum space then 
yields the following equations for the two-particle propagator 
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and the vertex function 

nkkdq, to) =2/ J' -^Gk +q (w + w')G_k(-<y')rkk'(<y. &/), 


(23) 


Tkk'(w, wO -<5kk' + ^ 2/ J' -^-Gk,+q (to + wOG-k^-Wl) 


X ^kk, ( ww l)rk,k'(w, Ctq), 


(24) 


with the basic one-magnon propagator up to 1 /S order defined 
as 


G- l (k,tj) = cj-oj k + iO + . (25) 

The factor of 2 in Eq.(23) and Eq.(24) stem from the two sets 
of contributions differing by the interchange of dummy mo¬ 
menta k(ki) + q and -k(kj) according to the Wick’s theorem. 
The lowest order two-particle irreducible interaction vertex 
Wir, shown in Fig. 5(c), reads as 


40]. An inspection of the interaction vertices for the TLAF 
reveals that r y^’ b) can not be separated into finite sum of prod¬ 
ucts of the triangular-lattice harmonics, thus Eq.(14) can not 
be algebraically solved in terms of a finite number of scat¬ 
tering channels. However, a numerical solution can be per¬ 
formed on finite lattices by summing over A points of k in 
the 1 st BZ, leading to a A x A system for the linear solver. 
We adopt this numerical approach to compute the interacting 
intensity plots. 

We begin with substituting (23) and (24) into (14), 

Xi = Y, /V/k/V/k' [dkk-1 Ik + I Ik ^ Vkk,n kl k-], (32) 

kk' k, 

where Hk = 2[o> - Sk+q - fc'k + /O ’]^ 1 represents the renor- 
malizated two-magnon propagator in the absence of vertex 
correction. The BZ on finite lattices can be divided into 
VA x vA meshes with the replacement of the continuous mo¬ 
menta (k,k',ki) to discrete variables (m,n, 1). The elements 
for the bimagnon susceptibility matrix are given by 


'Vm = 'Vq + + 'Vf + "Vf + 'Vf, (26) 

where the frequency-independent four-point vertex r V\ com¬ 
ing from the quartic Hamiltonian has the form 

W 4 = V c (ki + q, -ki; k + q, -k), (27) 

and the other four vertices q/ a_,/) j n the same 1/5 order 
which are assembled from two three-point vertices and one 
frequency-dependent propagator can be written as 

<Vf = 7^772 t y «( k i + q, k - k!; k + q)G 0 (k - kj, to' - to{) 
( 2 !) 

x K( _ k,k - k i; -kO], (28) 

*Vf = -^[V fl *(k + q, ki - k; k! + q)G 0 (k 1 - k, oq - to') 
(2!) 

xVa(-k|, k, - k; -k)], (29) 

^ + q.-kt;q)Go(q,w) 

(2!r 

xV*(k + q, -k; q)], (30) 

= 7 ^ 7 y[^(k 1 + q, -k l ,q)G 0 (-q, -to) 

(3!r 

xv;(k + q,-k,q)], (31) 

where we have retained only the bare propagator Go for each 
intermediate line in r V ( " ‘ l] in the spirit of 1/5 expansion. 
We further assume that two on-shell magnons are created 
and annihilated in the repeated ladder scattering process with 
to' ~ and <o t « -e^ 0) [51, 52], This approximation is 
best for sharp spectral peaks of the two main magnons in the 
scattering process where all the lowest order irreducible ver¬ 
tices are not explicitly frequency dependent. Based on the 
above simplifications, we now derive the final solution of the 
interacting RIXS intensity from the ladder approximation BS 
equation. 

An approach to solving the coupled BS equations is to 
decompose the irreducible vertices into lattice harmonics as 
demonstrated for the case of collinear antiferromagnet [38, 


X inn — .17 in A'7 n | () mnn m + n m Yj vunin]. (33) 

i 

We then obtain the eigenvalue equation for these discrete mo¬ 
menta 

^mn = Q mn z rmltfln, (34) 

1 

where the new functions are defined as 

^mn @mn ^mnnm^fn? l~ml — nm^ml■ (35) 

A direct solution to (34) gives the final form of the x matrix 
as 

x = m- fr 1 #, (36) 

where all the matrices in Eq. (36) have Ax A dimensions with 
the matrix elements explicitly defined as 

lmn — ^rnni^nm — 4 mn '17 nl , (37) 

f mn = n m y mn , Q mn = (5 mn n m M n . (38) 

The interacting pure bimagnon RIXS susceptibility can then 
be computed as 

Af 2 (q, to) = (39) 

m,n 

In Fig. 6 we plot the results for interacting bimagnon RIXS 
intensity. We choose two special BZ momenta values, M and 
Y, to illustrate our findings. In Fig. 6 (a) and Fig. 6 (d) we 
show the progression of the indirect RIXS spectra shape as 
the spin 5 value is changed from the classical case 5 —> oo 
(top), to 5 =3/2 (middle), to the maximal quantum case of 
5 = 1/2 (bottom). While the classical RIXS spectra from 
both momenta contain peaks due to the presence of Van Hove 
singularities, introduction of quantum fluctuations cause some 
of these spurious peaks to disappear. But observe that in the S 
= 1/2 case the spectra shape is strikingly different. In the ab¬ 
sence of anisotropy at the q = M (roton transfer momentum), 
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u/3 JS u/3JS u/3JS ui/3JS 


FIG. 6. (Color online) Interacting bimagnon RIXS intensity for transformed momenta q = M (a-c) and q = Y (d-f). In (a) and (d) the evolution 
of the interacting bimagnon intensity profile with increasing spin value 5 for the isotropic model is shown. In (a) and (d): 5 —> oo (top frame), 
5 = 3/2 (middle frame), and 5 = 1/2 (bottom frame). In (b) and (e) comparison of the contribution from the direct CV 4 ), indirect ( < V 3 ) and 
full (<V 3 + <V 4 ) vertices correction to the RIXS spectrum for 5 = 1/2 and a = 1 are displayed. In (c) and (f) the effects of easy-plane anisotropy 
on the splitting feature of the interacting bimagnon spectrum are shown. 


we observe a single peak at an energy of 4 JS . However, at the 
q —Y point there is a multipeak structure, see Fig. 6 (d). Now 
comparing with the magnon decay intensity map. Fig. 1(b), 
it is evident that there is a direct correlation between the sta¬ 
bility of the spin wave modes and the appearance of a single 
or multipeak structure. The above mentioned comparison is 
not restricted to these two choosen points. Comparision of 
the RIXS spectra generated from other special high symmetry 
momentum transfer also have the same features, see Fig. 7. 
Based on these observations we propose that RIXS can be 
used as a probe to detect the presence of the roton mode in 
a TLAF. Furthermore, to provide a comprehensive picture of 
the effects of geometrical frustration and anisotropy we intro¬ 
duce a small anisotropy a = 0.9 in the system. From Fig. 6 (c) 
and Fig. 6 (f) it is clear that inclusion of anisotropy causes fur¬ 
ther peak splitting. Thus a proper explanation of the RIXS 
spectra features in a TLAF involves analyzing both the effects 
of magnon damping and anisotropy. 

The 'V 4 and 'V 3 vertices play an important role in the gen¬ 
eration of the RIXS bimagnon spectra. Especially at the roton 
point it is worth noting that including only the direct colli¬ 
sion vertex does not renormalize the single-peak stmcture in 
the extreme quantum condition with S = 1/2, Fig. 6 (b). The 
major contribution to the interacting RIXS spectra originates 
from the indirect vertices arising from the three-magnon inter¬ 
action terms. This indicates that renormalization of the spectra 
is due to the indirect 'V 3 vertices which involve virtual decay 
and recombination of the two main magnons in the scattering 
process. This is different from the q = Y point where both the 
<V 4 and 'V 3 vertices contribute, as seen in Fig. 6 (e). 

Before we end this section it is important to point out an im¬ 
portant difference between a Raman scattering calculation and 
RIXS. In the case of RIXS, the contributions of diagrams 'V ( 3 C> 
and yf vanish identically when the transferred momenta be¬ 


(a) (b) 



FIG. 7. (Color online) Full indirect RIXS spectra U + / 3 of the 5 = 
1/2 XXZ triangular-lattice antiferromagnet with a = 1; (a) and (c) 
and a = 0.9; (b) and (d), along the high symmetry path T —> M and 
T —> K in the BZ. 

longs to the T — M path or related symmetrical lines in the BZ 
in accordance with the magnetic Raman scattering study for 
which q = 0 [51]. To demonstrate this fact we consider the 
contributions of the two-particle irreducible vertices 'V 3 and 
which are already in separated forms as functions of k 
and k'. The corresponding reducible vertex function Eq. (24) 
with respect to these diagrams can be directly obtained as 

r kk ' =<W + V* ( 6 ) (k)/(q, tu), (40) 

where /( q, at) is a function of to and q only. Barring the nonin- 
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FIG. 8. (Color online) Contour plot of the bimagnon velocity in the 
first BZ of TLAF for transferred momenta q = M (first row) and 
q = Y (second row) with a = 1 and a = 0.9. First column is the bi¬ 
magnon velocity in the absence of anisotropy. Second column is for 
the anisotropic (a = 0.9) case. Lowest velocity contours indicated 
by blue (grey) bands. Relative higher velocities indicated by black 
regions. 

teracting contributions the vertex correction to the RIXS sus¬ 
ceptibility is given by 

Xv(*F oj) = f( q, oj) ^ M k ' ^ M k n k R* (fc) (k) 

k' k 

= const x ^ M k U k V* (b) (k). (41) 

k 

In the above both M k and n k are even functions of k, while the 
function V* (bj (k) = V* Ah) {k + q, -k, q) are odd functions with 
respect to k when momentum q is along the FM line (e.g. q x = 
0). Thus, by virtue of the C 3 symmetry of hexagonal lattices 
we can conclude that the total contributions of diagrams 'V 1 ,' 1 
and 'V'-f 1 vanish identically when transformed momenta q are 
located in the lines from the center of the BZ to the middle 
of the BZ boundary. This implies that the source processes 
are prohibited in the repeated ladder scattering events when 
transferred momenta are along these symmetrical paths. 


C. Total RIXS intensity 

Using Eq. (13) we compute the full indirect RIXS spectra. 
In Fig. 7 we display the RIXS line plots along the F —> M 
path and along the T —> K path, respectively. The features ob¬ 
served are reminiscent of those discussed for the noninteract¬ 
ing trimagnon spectra and the full interacting bimagnon spec¬ 
tra. As noted earlier, we find that at the T point the spec¬ 
tra originates purely from the trimagnon contribution, irre¬ 
spective of the presence or absence of anisotropy. However, 
inclusion of anisotropy causes a downshift of the bimagnon 


contribution and an upward shift of the trimagnon spectra. 
Anisotropy gives rise to further splitting in the bimagnon case, 
however, the trimagnon spectra is not affected. The occurence 
of peak splitting observed in the RIXS spectrum can be pre¬ 
dicted by observing the bimagnon velocity plot. In a previ¬ 
ous publication on the square lattice Heisenberg magnet [40] 
we had highlighted the connection between bimagnon veloc¬ 
ity and the appearance of multipeak structure in the RIXS 
spectra. Interestingly enough, even within the context of a 
TLAF this relationship persists. To demonstrate this correla¬ 
tion, in Fig. 8 , we show the bimagnon velocity intensity plot 
in both the presence and absence of anisotropy for the q = M 
and the q - Y point. The black regions represent the high¬ 
est moving bimagnon velocities which clearly disappear with 
the inclusion of anistropy. As more puddles of slow moving 
bimagnon velocity appears, so does the appearance of a mul¬ 
tipeak structure as shown in Fig. 7. At the q = M point the 
single roton peak melts away with increasing anistropy which 
comes along with low bimagnon velocity. A similar effect is 
observed at the q = Y point, where with increasing anisotropy 
there are greater pockets of slow moving bimagnon. Hence, 
with anisotropy the peak splits further at Y point. Finally, in 
Fig. 9 we present the expected constant-ru scans of the to¬ 
tal interacting RIXS intensity for four selected energies from 
low energy to a high energy. One of the advantages of these 
constant-energy scans, which is reminiscent of the INS ex¬ 
periments, is that prominent peak structures are easy to dis¬ 
tinguish [61]. For the isotropic model the interacting RIXS 
intensity is strongly peaked on the corners of the hexagonal 
BZ at low and high energies, while these peaks disperse along 
the edges of the BZ at intermediate energies. However the 
presence of XXZ anisotropic strongly reduces these prominent 
features, which is in qualitative agreement with the noninter¬ 
acting total spectral weight shown in Fig. 2. 


V. CONCLUSION 

At present, there exists no theoretical guidance for experi¬ 
mentalists on how to analyze and interpret the RIXS spectra of 
an ordered phase in a geometrically frustrated quantum mag¬ 
net. Although a proposal for detecting spin-chirality terms in 
triangular lattice Mott insulators via RIXS has been put for¬ 
ward [ 66 ], there has been no analysis on the effect of geo¬ 
metrical frustration and anisotropy on the indirect RIXS spec¬ 
tra. In this paper, using a l/S expansion spin-wave theory 
involving Bethe-Salpeter corrections we investigate the key 
signatures of noncollinear ground state ordering in the indi¬ 
rect RIXS spectrum of a TLAF. We conclude that in the ab¬ 
sence of anisotropy the root cause of the multipeak structure 
is magnon decay. This mechanism is different from that of 
a square lattice where strong frustrating further neighbor in¬ 
teractions and anisotropy are required to cause peak splitting 
(instability). In the introduction we had put forward a cou¬ 
ple of questions- (a) How does the presence of an intrinsic 
damping affect the indirect K-edge RIXS spectra? and (b) 
What role does the interplay between geometrical frustration 
and spin anisotropy have on the RIXS spectra? Based on our 
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uj/3JS = 1.5 



uj/3JS = 2.0 


FIG. 9. (Color online) Intensity plots of the constant-w scans of the full RIXS intensity /(q, ai) scaled by (3 JS ) 2 for the 5 = 1/2 XXZ triangular 
lattice antiferromagnet in the q plane with a = 1 (upper) and a = 0.9 (lower) at four representative energies. The prominent peaks on the 
comers are strongly reduced in the presence of anisotropy in qualitative agreement with the noninteracting total spectral weight. 


calculations, we conclude that magnon damping does affect 
the spectra, causing the RIXS peak to be either stable (no 
splitting) or unstable (splitting leading to multipeak) in the 
absence or presence of damping, respectively. Geometrical 
frustration introduces noncollinear ordering which introduces 
magnon damping. The stability or instability of the ensuing 
magnon mode then dictates the appearance of a single or mul¬ 
tipeak structure. By comparing the K-edge RIXS intensity of 
the square lattice case, to that of the TLAF, we find that the 
RIXS intensity does not vanish at the F point and at the antifer¬ 
romagnetic wavevector. At the I point, the bimagnon inten¬ 
sity is zero and the single peak spectra results purely from the 
trimagnon contribution, approximately at energy scale of 6 JS 
corresponding to the three magnon energy. This provides ex¬ 
perimentalists with a means to detect purely trimagnon RIXS 
spectra at the K-edge. Our proposed scheme of detecting tri¬ 
magnon excitations is different from that put forward in the 
paper by Ament and Brink [44], since we are not consider¬ 
ing the L-edge. The single roton peak occurs at an energy 
of 4JS and can be used as an experimental signature to de¬ 
tect roton modes in a TLAF. In conclusion, our theoretical in¬ 
vestigation demonstrates that RIXS has the potential to probe 
and provide a comprehensive characterization of the micro¬ 
scopic properties of bimagnon and trimagnon excitations in 
the TLAF across the entire BZ, which is beyond the capabili¬ 
ties of traditional low-energy optical techniques [49-52], 
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Appendix A: Derivation of interacting spin-wave theory 


We utilize the Holstein-Primakoff transformation to 
bosonize the local rotating Hamiltonian (1) 

SJ = S - a-a,, SJ = a f ^2 S -ate?,, S, + = (SJ)\ (Al) 

with subsequent expansion of square root to first order in 
a\cijl2S. This is followed by a Fourier transformation. The 
Fourier transformed Hamiltonian takes the form 

<H = H 0 + H 2 + H 3 + H 4 + 0(S “'). (A2) 


The first term corresponds to the classical energy and the 
quadratic Hamiltonian reads 

Ho = Y + ^B k (a[fli k + a- k fl k )], 

k 

A k = 3 JS [1 + (a - ^)y k ], B k = -3 JS(a + '-)y k , (A3) 
with the structure factor y k defined as 


h y k * 

7k = cos k x + 2 cos — cos —k y ). 


(A4) 


We then diagonalize the harmonic part /7 2 by the Bogoliubov 
transformation 


fl k = u k b k + v k b_ k , (A5) 

with the parameters t/ k and i; k defined as 

, 2 A k + e k B k 

u k’ v k = —3 - ,u k v k = ~—, (A6) 

2 Ek 2e k 
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and the linear spin-wave theory dispersion given by 

£ k = 7 A k ~ B l = 3^'e k , (A7) 

where we have defined the dimensionless energy 

£ k = y/(l - yk)(l + 2y k ). (A8) 

Performing the Bogoliubov transformations in the cubic in¬ 
teraction term H 3 we obtain 

H 3=^ X K,(ki,k 2 ;k 3 )(^ i ^/7 k3 +H.c.) 

k] +k 2 =k 3 

+ h 2 n(k 1 .k 2 ,k 3 )(^ i ^ 2 ^ 3 +H.c.)(A9) 

k]+k 2 +k 3 =0 


t 2 a + 1 

5 = ^ = — 4 — Cu ( A16 ) 

with the two-dimensional integrals 

(y k V 

Cl (l = 0,1,2). (A17) 

k ^ 

After the mean-field decoupling, the quartic part is decom¬ 
posed as 

H 4 = 5H 0 + 6H 2 + Ha. (A18) 

The first term 6 Hq is the correction to the ground-state energy 
and the quadratic parts reads 


The explicit forms for the three-boson interaction vertices are 


Va(l,2;3) =37ty — \y\(u\ + i>i )(u 2 u 3 + v 2 v 3 ) + y 2 (u 2 + v 2 ) 
x(u x u 3 + uiu 3 ) - y 3 (u 3 + v 3 )(u x v 2 + i>i« 2 )],(A10) 
3S", 


14(1,2,3) =3/t yj — [yi(ni + vi )(u 2 v 3 + v 2 u 3 ) + y 2 (u 2 + v 2 ) 
x(m iw 3 + yi« 3 ) + y 3 (u 3 + w 3 )(miw 2 + dim 2 )],(A1 1) 


where u,, o,- are Bogoliubov parameters and the function y k is 
defined as 


yk =!(sin£ x -2siny cos—(A12) 


The three-boson vertex V a and 14 in // 3 describes interaction 
between one- and two-magnon states and are called the decay 
and the source vertex, respectively. 

To derive the explicit forms of the quartic interaction term 
Ha , it is convenient to introduce the following Hartree-Fock 
averages 



1 

2 a- 1 1 


n = (a cii) — 

2 C0 + 

4 Cl _ 2’ 

(A13) 


1 

2 a- 1 


m = {a a,j) — 

2 C1 + 

4 

(A14) 


2a + 1 



A = ( didj ) = 

4 

c 2 . 

(A15) 


5H 2 = V [&4 k a£a k + ] -6B k (alci r k + a_ k a k )], (A19) 

k 

with 

&4 k =^[(1 + 2a)A + (1 - 2 a)m - 2n\ 

3 

+ -[(1 + 2a)6 + 2(1 - 2 a)n - 4/n]y k , 

AB k =-[(1 + 2 a)m + (1 - 2a)A] 

+ ^[(1 - 2a)<5 + 2(+2a)«-4A]y k . (A20) 

We then obtain the Hartree-Fock correction to the harmonic 
spin-wave spectrum 


fe k = (« k + u k )&4 k + 2 M k y k hB k . (A21) 

The normal-ordered term H 4 describes the multi-particle in¬ 
teractions. Here we only display the explicit expression for 
the lowest order irreducible two-particle scattering amplitude 
which is relevant for the our calculations as 

H\ P = 2 K.(k,,k 2 ;k 3 ,k 4 )^ i y/7 k3 /7 kj ,(A22) 

k [ +k 2 =k 3 +k 4 


with the vertex function 


3 JS i 

14(1,2; 3,4) =Y^r{(2a + l)(y! +y 2 + r 4 )(u 1 u 2 u 3 V4 + v 1 v 2 v 3 u 4 ) + (2a + l)(yi + y 2 + y 3 )(u x u 2 v 3 u 4 + v x v 2 u 3 v 4 ) 

+(2 a + l)(y 2 + y 3 + y 4 )(uiv 2 u 3 u 4 + viu 2 v 3 v 4 ) + (2a + l)(yi + y 3 + y 4 )(iti 1^0314 + v\u 2 u 3 u 4 ) 

—[ 2 (yi_ 3 + y 2 - 3 + yi _4 + y 2 - 4 ) + (2a - l)(yi + y 2 + y 3 + y 4 )](u x u 2 u 3 u 4 + v x v 2 v 3 v 4 ) 

— [ 2 (yi +2 + y 3+4 + yi _ 3 + y 2 _ 4 ) + (2a - lXyr + y 2 + y 3 + y 4 )](u\v 2 u 3 v 4 + v x u 2 v 3 u 4 ) 

-[ 2 (yi +2 + y 3+4 + yi - 4 + yi- 3 ) + (2a - l)(yi + y 2 + y 3 + y 4 )](uiv 2 v 3 u 4 + y 1 M 2 n 3 y 4 )}. (A23) 


By collecting all these terms together, we finally obtain the Appendix B: Exact versus numerical solution to the BS equation 
effective interacting spin-wave Hamiltonian Eq. (2). 

To test the validity of our numerical method on finite lat¬ 
tices (N — 69 x 69), we adopt the exact solution approach to 
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TABLE 1. Definition of the channels u n (k) 


n 

w„(k) 

n 

f«(k) 

i 

M k+q w k 

15 

^k+q^k COS k x COS ky 

2 

fk+qt'k 

16 

Vk+qUk COS k x COS -y- ky 

3 

^k+q^k 

17 

Uk+qUk sinfe v sin -yk y 

4 

^k+q w k 

18 

fk+q^k sin k x sin -yk y 

5 

Uk+qMk COS k x 

19 

iik+qVk sin& Y sin -yk y 

6 

f'k+ql'k COS k x 

20 

Vk+qUk sin ^ sin -yk y 

7 

U k+qUk COS k x 

21 

Hk+q/<k cos k x sin -yk y 

8 

^k+q^k COS k x 

22 

t’k+q^k cos k x sin y k tJ 

9 

u k+q u k sin k x 

24 

Uk+qVk cos k x sin y k tJ 

10 

^k+q ^k $in k x 

24 

Vk+qUk cos k x sin -yk y 

11 

^k+q^k SHI k x 

25 

Uk+qUk sin£ r COS yk y 

12 

^k+q ^k sin k x 

26 

t’k+q^k sin k x COS yky 

13 

Uk+qUk COS k x COS ~yky 

27 

Uk+qVk sin& Y cos y k tJ 

14 

^k+q^k COS k x COS ky 

28 

^k+qi'k sin k x cos ky 


solving a BS equation outlined in Appendix B of our publica¬ 
tion Ref. [40]. We obtain a separated form for the four-point 
vertex 'V 4 for the Heisenberg model (a = 1) on triangular 
lattice which has the following expression 


28 

*V 4 (ki + q, —ki; k + q, -k) = ^ y m (k)f mn u n (ki). (Bl) 

m,n=l 



The channels u„(k) are defined in Table I with the matrix ele¬ 
ments of f denoted by 


3 JS § i f \ 

165 f §2 ) 


(B2) 


FIG. 10. (Color online) Interacting bimagnon RIXS intensity renor¬ 
malized only by the direct interaction vertex (TV) based on the exact 
(N —> oo) versus numerical (TV = 69 x 69) solution of the Bethe- 
Salpeter equation for a S = 1/2 isotropic triangular lattice antiferro- 

magnet at various momenta in the BZ. where the blocks ar£ given by 
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'-§ -|/i o o o o o o §4 o o o |v ' 

-§ 0 0 -§* 0 0 0 fyu 0 0 0 fv 0 

-§0 0 0 0 -fv 0 0 0 -fju 0 0 0 

-§0 0 0 0 -fv 0 0 0 -fyu 0 0 

-f -§T 0 0 0 -fv 0 -fyu 0 0 

-f 0 0 -fv 0 0 0 -§// 0 

y 0 0 0 -\x 0 0 0 

fv 0 0 0 -lx 0 0 

-f -§T 0 0 0 -\x 

-f 0 0 -lx 0 

\tp 0 0 0 



In the above we have introduced the following notations: 

C° = cos q x , c\ = cos y, C 2 3 4 5 6 7 8 = cos ^-q y , (B3) 

5° = cos q x , S\ = cos y, 5^ = cos ^-q y , (B4) 

A = C l q C 2 q , n = C\S 2 q , v = s\c\,x = S\S 2 q , (B5) 
0 = C° + 1, t? = C°-l, (B6) 

0 = C‘C* + 1, y = CjC*-l. (B7) 


Only the upper right parts of S i and S 2 are shown since the 
matrices are symmetrical. 


In Fig. 10 we show the results of our computation for 
S = 1/2 and a — 1. The numerical solution of Eq. (39) by 
considering only the 'V 4 vertex is performed with 69 x 69 lat¬ 
tice sites, while the integrals arising in the exact solution are 
solved on a mesh of size 252 x 252 and then extrapolated to 
IV —> 00 [40]. 


[1] P. W. Anderson, “Resonating valence bonds: A new kind of 
insulator?” Mater. Res. Bull. 8, 153 (1973). 

[2] Th. Jolicoeur and J. C. Le Guillou, “Spin-wave results for the 
triangular heisenberg antiferromagnet,” Phys. Rev. B 40, 2727- 
2729 (1989). 

[3] Satoru J. Miyake, “Spin-wave results for the staggered magne¬ 
tization of triangular heisenberg antiferromagnet,” J. Phys. Soc. 
Jpn. 61, 983-988 (1992). 

[4] A V Chubukov, S Sachdev, and T Senthil, “Large-s expansion 
for quantum antiferromagnets on a triangular lattice,” J. Phys.: 
Condens. Matter 6, 8891 (1994). 

[5] Luca Capriotti, Adolfo E. Trumper, and Sandro Sorella, “Long- 
range neel order in the triangular heisenberg model,” Phys. Rev. 
Lett. 82, 3899-3902 (1999). 

[6] Weihong Zheng, John O. Fjasrestad, Rajiv R. P. Singh, Ross H. 
McKenzie, and Radu Coidea, “Excitation spectra of the spin- 3 
triangular-lattice heisenberg antiferromagnet,” Phys. Rev. B 74, 
224420 (2006). 

[7] Steven R. White and A. L. Chernyshev, “Neel order in square 
and triangular lattice heisenberg models,” Phys. Rev. Lett. 99, 
127004 (2007). 

[8] P. H. Y. Li, R. F. Bishop, and C. E. Campbell, “Quasiclassical 
magnetic order and its loss in a spin- / heisenberg antiferromag¬ 
net on a triangular lattice with competing bonds,” Phys. Rev. B 
91,014426(2015). 


[9] David A. Huse and Veit Elser, “Simple variational wave func¬ 
tions for two-dimensional heisenberg spin-'/i antiferromag¬ 
nets,” Phys. Rev. Lett. 60, 2531-2534 (1988). 

[10] B. Bemu, P. Lecheminant, C. Lhuillier, and L. Pierre, “Exact 
spectra, spin susceptibilities, and order parameter of the quan¬ 
tum heisenberg antiferromagnet on the triangular lattice,” Phys. 
Rev. B 50, 10048-10062 (1994). 

[11] Rajiv R. P. Singh and David A. Huse, “Three-sublattice order 
in triangular- and kagome-lattice spin-half antiferromagnets,” 
Phys. Rev. Lett. 68, 1766-1769 (1992). 

[12] Yutaka Shirata, Hidekazu Tanaka, Akira Matsuo, and Koichi 
Kindo, “Experimental realization of a spin-1/2 triangular- 
lattice heisenberg antiferromagnet,” Phys. Rev. Lett. 108, 
057205 (2012). 

[13] G. Koutroulakis, T. Zhou, Y. Kamiya, J. D. Thompson, H. D. 
Zhou, C. D. Batista, and S. E. Brown, “Quantum phase diagram 
of the s = \ triangular-lattice antiferromagnet ba 3 Cosb 209 ,” 
Phys. Rev. B 91, 024410 (2015). 

[14] Takuya Susuki, Nobuyuki Kurita, Takuya Tanaka, Hi- 
royuki Nojiri, Akira Matsuo, Koichi Kindo, and Hidekazu 
Tanaka, “Magnetization process and collective excitations 
in the i=l/2 triangular-lattice heisenberg antiferromagnet 
ba 3 cosb 2 09 ,” Phys. Rev. Lett. 110, 267201 (2013). 

[15] T. Ono, H. Tanaka, H. Aruga Katori, F. Ishikawa, H. Mitamura, 
and T. Goto, “Magnetization plateau in the frustrated quantum 






14 


spin system cs 2 cubr 4 Phys. Rev. B 67, 104431 (2003). 

[16] Hiroaki Kadowaki, Koji Ubukoshi, Kinshiro Hirakawa, Jos L. 
Martnez, and Gen Shirane, “Experimental study of new type 
phase transition in triangular lattice antiferromagnet vcl2,” J. 
Phys. Soc. Jpn 56, 4027-4039 (1987). 

[17] R. Ishii, S. Tanaka, K. Onuma, Y. Nambu, M. Tokunaga, 
T. Sakakibara, N. Kawashima, Y. Maeno, C. Broholm, D. P. 
Gautreaux, J. Y. Chan, and S. Nakatsuji, “Successive phase 
transitions and phase diagrams for the quasi-two-dimensional 
easy-axis triangular antiferromagnet rb 4 mn(moo 4 ) 3,” Eur. 
Phys. Lett. 94. 17001 (2011). 

[18] M. Poienar, F. Damay, C. Martin, J. Robert, and S. Petit, “Spin 
dynamics in the geometrically frustrated multiferroic cucro 2 ,” 
Phys. Rev. B 81, 104411 (2010). 

[19] S. Toth, B. Lake, S. A. J. Kimber, O. Pieper, M. Reehuis, 
A. T. M. N. Islam, O. Zaharko, C. Ritter, A. H. Hill, H. Ryll, 
K. Kiefer, D. N. Argyriou, and A. J. Williams, “120° heli¬ 
cal magnetic order in the distorted triangular antiferromagnet 
a-ca 2 cubr 4 Phys. Rev. B 84, 054452 (2011). 

[20] S. Toth, B. Lake, K. Hradil, T. Guidi, K. C. Rule, M. B. Stone, 
and A. T. M. N. Islam, “Magnetic soft modes in the distorted 
triangular antiferromagnet a-ca 2 cubr 4 ,” Phys. Rev. Lett. 109, 
127203 (2012). 

[21] R. Coidea, D. A. Tennant, A. M. Tsvelik, and Z. Tylczynski, 
“Experimental realization of a 2d fractional quantum spin liq¬ 
uid,” Phys. Rev. Lett. 86 , 1335-1338 (2001). 

[22] Ru Chen, Hyejin Ju, Hong-Chen Jiang, Oleg A. Starykh, and 
Leon Balents, “Ground states of spin - 3 triangular antiferromag- 
nets in a magnetic field,” Phys. Rev. B 87, 165123 (2013). 

[23] Burkhard Schmidt and Peter Thalmeier, “Quantum fluctuations 
in anisotropic triangular lattices with ferromagnetic and antifer¬ 
romagnetic exchange,” Phys. Rev. B 89, 184402 (2014). 

[24] Philipp Hauke, Tommaso Roscilde, Valentin Murg, J Ignacio 
Cirac, and Roman Schmied, “Modified spin-wave theory with 
ordering vector optimization: spatially anisotropic triangular 
lattice and j 1 j 2 j 3 model with heisenberg interactions,” New 
Journal of Physics 13, 075017 (2011). 

[25] Masanori Kohno, Oleg A. Starykh, and Leon Balents, “Spinons 
and triplons in spatially anisotropic frustrated antiferromag- 
nets,” Nature Phys. 3, 790 (2007). 

[26] M. Swanson, J. T. Haraldsen, and R. S. Fishman, “Critical 
anisotropies of a geometrically frustrated triangular-lattice an¬ 
tiferromagnet,” Phys. Rev. B 79, 184413 (2009). 

[27] Randy S. Fishman and Satoshi Okamoto, “Noncollinear mag¬ 
netic phases of a triangular-lattice antiferromagnet and of doped 
cufeo 2 ,” Phys. Rev. B 81, 020402 (2010). 

[28] E. A. Ghioldi, A. Mezio, L. O. Manuel, R. R. P. Singh, J. Oit- 
maa, and A. E. Trumper, “Magnons and excitation contin¬ 
uum in xxz triangular antiferromagnetic model: Application to 
ba 3 cosb 2 o 9 ,” Phys. Rev. B 91, 134423 (2015). 

[29] Nobuo Suzuki, Fumitaka Matsubara, Sumiyoshi Fujiki, and 
Takayuki Shirakura, “Absence of classical long-range order in 
an s = 5 heisenberg antiferromagnet on a triangular lattice,” 
Phys. Rev. B 90, 184414 (2014). 

[30] Andreas Weichselbaum and Steven R. White, “Incommensu¬ 
rate correlations in the anisotropic triangular heisenberg lat¬ 
tice,” Phys. Rev. B 84, 245130 (2011). 

[31] Philipp Hauke, “Quantum disorder in the spatially completely 
anisotropic triangular lattice,” Phys. Rev. B 87, 014415 (2013). 

[32] R. Coidea, S. M. Hayden, G. Aeppli, T. G. Perring, C. D. Frost, 
T. E. Mason, S.-W. Cheong, and Z. Fisk, “Spin waves and elec¬ 
tronic interactions in la 2 cuo 4 ,” Phys. Rev. Lett. 86 , 5377-5380 
( 2001 ). 


[33] H. M. Rpnnow, D. F. McMorrow, R. Coidea, A. Harrison, I. D. 
Youngson, T. G. Perring, G. Aeppli, O. Syljuasen, K. Lefmann, 
and C. Rischel, “Spin dynamics of the 2d spin I quantum anti¬ 
ferromagnet copper deuteroformate tetradeuterate (cftd),” Phys. 
Rev. Lett. 87, 037202 (2001). 

[34] Luuk J. P. Ament, Michel van Veenendaal, Thomas P. Dev- 
ereaux, John P. Hill, and Jeroen van den Brink, “Resonant in¬ 
elastic x-ray scattering studies of elementary excitations,” Rev. 
Mod. Phys. 83, 705-767 (2011). 

[35] J. P. Hill, G. Blumberg, Young-June Kim. D. S. Ellis, S. Waki- 
moto, R. J. Birgeneau, Seiki Korniya, Yoichi Ando, B. Liang, 
R. L. Greene, D. Casa, and T. Gog, “Observation of a 500 mev 
collective mode in la 2 _ r sr v cuo 4 and nd 2 cuo 4 using resonant in¬ 
elastic x-ray scattering,” Phys. Rev. Lett. 100, 097001 (2008). 

[36] D. S. Ellis, Jungho Kim, J. P. Hill, S. Wakimoto, R. J. Bir¬ 
geneau, Y. Shvyd’ko, D. Casa, T. Gog, K. Ishii, K. Ikeuchi, 
A. Paramekanti, and Young-June Kim, “Magnetic nature of the 
500 mev peak in la 2 _ x sr A cuo 4 observed with resonant inelas¬ 
tic x-ray scattering at the Cufc-edge,” Phys. Rev. B 81, 085124 
( 2010 ). 

[37] J. van den Brink, “The theory of indirect resonant inelastic x- 
ray scattering on magnons,” Europhys. Lett. 80, 47003 (2007). 

[38] Tatsuya Nagao and Jun-Ichi Igarashi, “Two-magnon excitations 
in resonant inelastic x-ray scattering from quantum heisenberg 
antiferromagnets,” Phys. Rev. B 75, 214414 (2007). 

[39] Filomena Forte, Luuk J. P. Ament, and Jeroen van den Brink, 
“Magnetic excitations in la 2 Cuo 4 probed by indirect resonant 
inelastic x-ray scattering,” Phys. Rev. B 77, 134428 (2008). 

[40] Cheng Luo, Trinanjan Datta, and Dao Xin Yao, “Spec¬ 
trum splitting of bimagnon excitations in a spatially frustrated 
heisenberg antiferromagnet revealed by resonant inelastic x-ray 
scattering,” Phys. Rev. B 89, 165103 (2014). 

[41] C. J. Jia, C.-C. Chen, A. P. Sorini, B. Moritz, and T. P. De- 
vereaux, “Uncovering selective excitations using the resonant 
profile of indirect inelastic x-ray scattering in correlated mate¬ 
rials: observing two-magnon scattering and relation to the dy¬ 
namical structure factor,” New Journal of Physics 14, 113038 
( 2012 ). 

[42] Luuk J. P. Ament, Giacomo Ghiringhelli, Marco Moretti Sala, 
Lucio Braicovich, and Jeroen van den Brink, “Theoretical 
demonstration of how the dispersion of magnetic excitations in 
cuprate compounds can be determined using resonant inelastic 
x-ray scattering,” Phys. Rev. Lett. 103, 117003 (2009). 

[43] M. W. Haverkort, “Theory of resonant inelastic x-ray scatter¬ 
ing by collective magnetic excitations,” Phys. Rev. Lett. 105, 
167404 (2010). 

[44] Luuk J. P. Ament and Jeroen van den Brink, “Strong 
three-magnon scattering in cuprates by resonant x-rays,” 
arXiv: 1002.3773. 

[45] M. E. Zhitomirsky and A. L. Chernyshev, “Colloquium: Spon¬ 
taneous magnon decays,” Rev. Mod. Phys. 85, 219-242 (2013). 

[46] A. L. Chernyshev and M. E. Zhitomirsky, “Magnon decay in 
noncollinear quantum antiferromagnets,” Phys. Rev. Lett. 97, 
207202 (2006). 

[47] J. van den Brink and M. van Veenendaal, “Correlation func¬ 
tions measured by indirect resonant inelastic x-ray scattering,” 
Europhys. Lett. 73, 121 (2006). 

[48] Luuk J. P. Ament, Filomena Forte, and Jeroen van den Brink, 
“Ultrashort lifetime expansion for indirect resonant inelastic x- 
ray scattering,” Phys. Rev. B 75, 115118 (2007). 

[49] Thomas P. Devereaux and Rudi Hackl, “Inelastic light scatter¬ 
ing from correlated electrons,” Rev. Mod. Phys. 79, 175-233 
(2007). 


15 


[50] F Vernay, T P Devereaux, and M J P Gingras, “Raman scat¬ 
tering for triangular lattices spin- 1/2 heisenberg antiferromag- 
nets,” J. Phys.: Condens. Matter 19, 145243 (2007). 

[51] Natalia Perkins and Wolfram Brenig, “Raman scattering in a 
heisenberg s = \ antiferromagnet on the triangular lattice,” 
Phys. Rev. B 77, 174412 (2008). 

[52] Natalia B. Perkins, Gia-Wei Chem, and Wolfram Brenig, “Ra¬ 
man scattering in a heisenberg s = 1 antiferromagnet on the 
anisotropic triangular lattice,” Phys. Rev. B 87, 174423 (2013). 

[53] Weihong Zheng, John O. Fjterestad, Rajiv R. P. Singh, Ross H. 
McKenzie, and Radu Coidea, “Anomalous excitation spectra 
of frustrated quantum antiferromagnets,” Phys. Rev. Lett. 96, 
057201 (2006). 

[54] Oleg A. Starykh, Andrey V. Chubukov, and Alexander G. 
Abanov, “Flat spin-wave dispersion in a triangular antiferro¬ 
magnet,” Phys. Rev. B 74, 180403 (2006). 

[55] A. L. Chernyshev and M. E. Zhitomirsky, “Spin waves in a tri¬ 
angular lattice antiferromagnet: Decays, spectrum renormaliza¬ 
tion, and singularities,” Phys. Rev. B 79, 144416 (2009). 

[56] H. D. Zhou, Cenke Xu, A. M. Hallas, H. J. Silverstein, C. R. 
Wiebe, I. Umegaki, J. Q. Yan, T. P. Murphy, J.-H. Park. Y. Qiu, 
J. R. D. Copley, J. S. Gardner, and Y. Takano, “Successive 
phase transitions and extended spin-excitation continuum in the 

triangular-lattice antiferromagnet basCosbjOg,” Phys. Rev. 
Lett. 109, 267206 (2012). 

[57] R. Coidea, D. A. Tennant, and Z. Tylczynski, “Extended scat¬ 
tering continua characteristic of spin fractionalization in the 
two-dimensional frustrated quantum magnet csiCucL observed 
by neutron scattering,” Phys. Rev. B 68 , 134424 (2003). 


[58] Joosung Oh, Manh Due Le, Jaehong Jeong, Jung-hyun Lee, 
Hyungje Woo, Wan-Young Song, T. G. Perring, W. J. L. Buy¬ 
ers, S.-W. Cheong, and Je-Geun Park, “Magnon breakdown in 
a two dimensional triangular lattice heisenberg antiferromagnet 
of multiferroic lunino^,” Phys. Rev. Lett. Ill, 257202 (2013). 

[59] M. Y. Veillette, A. J. A. James, and F. H. L. Essler, “Spin dy¬ 
namics of the quasi-two-dimensional spin-^ quantum magnet 
cs 2 cucl 4 ,” Phys. Rev. B 72, 134429 (2005). 

[60] Denis Dalidovich, Rastko Sknepnek, A. John Berlinsky, Jun- 
hua Zhang, and Catherine Kallin, “Spin structure factor of the 
frustrated quantum magnet CS 2 CUCI 4 ,” Phys. Rev. B 73, 184403 
(2006). 

[61] M. Mourigal, W. T. Fuhrman, A. L. Chernyshev, and M. E. 
Zhitomirsky, "Dynamical structure factor of the triangular- 
lattice antiferromagnet,” Phys. Rev. B 88 , 094407 (2013). 

[62] C. M. Canali and Mats Wallin, “Spin-spin correlation functions 
for the square-lattice heisenberg antiferromagnet at zero tem¬ 
perature,” Phys. Rev. B 48, 3264-3280 (1993). 

[63] J. Lorenzana, G. Seibold, and R. Coidea, “Sum rules and 
missing spectral weight in magnetic neutron scattering in the 
cuprates,” Phys. Rev. B 72, 224511 (2005). 

[64] R. W. Davies, S. R. Chinn, and H. J. Zeiger, “Spin-wave ap¬ 
proach to two-magnon raman scattering in a simple antiferro¬ 
magnet,” Phys. Rev. B 4, 992-1004 (1971). 

[65] C. M. Canali and S. M. Girvin, “Theory of raman scattering in 
layered cuprate materials,” Phys. Rev. B 45, 7127-7160 (1992). 

[ 66 ] Wing-Ho Ko and Patrick A. Lee, “Proposal for detecting spin- 
chirality terms in mott insulators via resonant inelastic x-ray 
scattering,” Phys. Rev. B 84, 125102 (2011). 


